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A.  INTRODUCTION 

Within  the  past  several  decades,  interest  in  the 
study  of  infra-red  and  Raman  spectra  has  erown  considerably. 
When  Rawlins  and  Taylor  published  their  book"*-  on  infra-red 
analysis  in  1928  ^  there  wasn't  a  single  book  available  con- 
cerned with  infra-red  study  of  molecular  structure.  Since 
then  a  number  of  books  have  appeared,  and  many  articles  have 
been  published  on  the  subject. 

There  are  several  reasons  for  the  growth  of  inter- 
est.     One  of  the  important  reasons  is  the  application  of 
infra-red  study  to  chemical  analysis.     It  is  possible  by  an  . 
analysis  of  the  infra-red  spectrum  of  a  substance  not  only 
to  identify  elements  but  atomic  and  molecular  croups  as  well. 
For  example,  an  ethyl  or  a  methyl  group  in  a  comoound  will 
always  exhibit  certain  characteristic  frequencies. 

Secondly,  the  study  of  molecular  spectra  has  proba- 
bly become  the  most  important  means  for  investigating  mole- 
cular structure.    Various  energy  levels  may  be  obtained  from 
a  study  of  the  spectra  which,  in  turn,  prive  detailed  infor- 
mation on  the  electronic  motion,  vibration,  and  rotation  of 
the  atoms  in  the  molecules.    A  theoretical  understanding  of 
valence  has  resulted  from  a  study  of  electronic  motion  while 
the  forces  acting  between  atoms  in  a  molecule  may  be  ob- 
tained from  the  vibrational  frequencies  and  accurate  infor- 


1.      Rawlins,  F.  I.  G. ,  and  Taylor,  A.M.,  Infra-Red  Analysis 

of  Molecular  Structure  (Cambridge  University  Press,  1929) 


ii 


mation  concerning  the  geometrical  arrangement  of  the  atoms 
may  be  obtained  from  a  study  of  the  rotational  spectra. 

Furthermore,  the  Raman  effect  was  discovered  in 
1928  which  also  has  become  important  in  the  field  of  analy- 
sis and  theoretical  studies.     It  is  sometimes  possible  ^o 
obtain  information  about  a  molecule  from  the  Raman  spectrum 
which  cannot  be  obtained  by  other  means. 

The  theoretical  aspects  of  the  study  of  molecular 
structure  made  considerable  strides  forward  when  group  theory 
was  introduced.    The  determination  of  the  structure  of  poly- 
atomic molecules  and  the  forces  acting  between  the  atoms  is 
an  extremely  complex  problem  but  a  great  deal  of  simplifi- 
cation is  introduced  by  the  application  of  grouo  theory  to 
the  symmetry  of  the  molecule. 

The  advance  of  theoretical  methods  always  requires 
an  improvement  of  experimental  techniques.     The  improved  ex- 
perimental methods  then  uncover  further  information  for  use 
in  theoretical  study.     As  a  consequence  of  improved  methods 
in  both  endeavors,  much  work  is  being  done  at  the  present 
time  in  molecular  studies.     It  is  with  the  theoretical 
aspects  of  the  problem  that  this  paper  is  concerned. 

The  infra-red  spectrum  is  divided  into  three 
general  regions  called  the  photoelectric,  near  infra-red 
and  far  infra-red  .     The  photoelectric  region  extends  from 

2.      Harrison,  George  R..  Lord,  Richard  C,  and  Loofbourow, 

John  R#,  Practical  Spectroscopy  (New  York:  Prentice-Kali 
Inc..  194F5 
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about  0.75  microns  to  about  2  microns  while  the  near  infra- 
red extends  from  about  2  microns  to  about  25  microns  and 
the  far  infra-red  from  about  25  microns  to  about  1  mm.  It 
is  divided  into  these  regions  because  of  the  experimental 
techniques  involved.     In  the  photoelectric  portion,  ordinary- 
prisms,  gratings  and  photographic  technique  are  utilized. 
However,  in  the  near  and  far  infra-red,  the  detection  of  the 
spectra  is  primarily  accomplished  b,r  utilizin-  the  heating 
effect  of  the  radiant  energy.     In  the  near  infra-red  prisms 
are  still  used,  but  they  are  made  of  rock  salt  or  potassium 
bromide,  while  in  the  far  infra-red.  gratings  are  used  prima- 
rily. 

These  divisions  in  a  general  way  agree  with  the 
division  of  the  molecular  energy  levels  into  electronic, 
vibrational  and  rotational  regions.     Infra-red  in  general 
corresponds  to  oscillations  of  particles  which  are  several 
hundred  times  slower  than  the  motion  of  electrons.  This 
rotion  is  attained  by  the  atoms  in  a  molecule. 

The  Raman  effect  is  a  scattering  process  which 
differs  from  the  Tyndall  ef fecit  in  that  with  the  Raman 
effect  there  is  a  change  of  frequency  associated  wi^h  the 
scattering.    Th<=>  incident  quanta  collide  with  a  vibrating 
atom  and  are  scattered  with  either  a  higher  or  lower  energy, 
depending  on  the  nature  of  the  collision,  and  hence  will 
have  a  resultant  higher  or  lower  frequency.     The  change  in 
frequency  is  characteristic  of  the  frequency  of  vibration 


of  the  atom  with  which  fhe  quantum  of  light  collides.  The 
Raman  lines  will  then  appear  in  a  position  displaced  from 
the  incident  spectral  line  by  an  amount  which  is  character- 
istic for  the  medium  through  which  the  light  has  been 
scattered.     The  lines  which  are  displaced  to  the  longer 
wave  length  side  are  called  Stokes  lines  and  those  displaced 
to  shorter  wave  lengths  are  called  anti-Stokes  after  a 
similar  phenomenon  which  appears  in  fluorescence  and  was 
discovered  b*r  Smokes.     In  fluorescence,  however,  there  is 
an  actual  absorption  of  the  light  quanta  which  after  a  short 
time  are  reemitted  with  a  change  of  wave  length.     Since  a 
relatively  small  amount  of  the  light  quanta  will  collide 
with  atoms,  the  intensity  of  the  Raman  lines  will  be  rather 
low  and  hence  the  experimental  difficulties  in  observing 
the  lines  are  correspondingly  greater. 

The  vibration  of  the  atoms  in  a  moleucle  give  rise 
+  o  the  major  portion  of  the  near  infra-red  and  to  Raman  lines 
A  knowledge  of  the  vibration  of  the  atoms  will  yield  infor- 
mation concerning  the  structure  of  the  molecule  and  par- 
ticularly the  forces  acting  between  the  atoms ' in  a  molecule. 
These  can  then  be  used  to  calculate  +he  energy  of  dissoci- 
ation and  other  thermodynamic  quantities  for  the  molecule, 
such  as,  heat  capacities,  etc.     It  is  with  this  facet  of  the 
problem,  that  is,  a  determination  of  the  normal  frequencies 
of  vibration  of  1,  3"  butadiene,  that  we  concern  ourselves 
primarily  in  this  paper. 
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B.-      THEORY  OF  MOLECULAR  VIBRATIONAL  SPECTRA 

1.  Molecular  Structure 

The  atoms  in  a  molecule  are  continually  in  motion 
vibrating  about  an  equilibrium  position  and  rotating  about 
their  center  of  mass  as  long  as  the  molecule  possesses  any 
thermal  energy.    The  equilibrium  configuration  is  best  de- 
termined from  a  study  of  the  rotational  spectra  or  electron 
diffraction  patterns.     The  vibrations  of  the  atoms  are  then 
considered  in  terms  of  the  displacements  from  these  equi- 
librium configurations. 

A  knowledge  of  the  number  and  size  of  the  atoms 
and  their  configuration  permits  a  utilization  of  their 
symmetry  properties  in  the  determination  of  the  frequencies 
of  vibration.     The  particular  configuration  which  the  atoms 
will  be  in  is  determined  by  the  forces  which  are  acting 
between  them.     These  forces  also  determine  the  normal 
frequencies  of  vibration. 

2.  Interatomic  Forces 

There  are  essentially  two  types  of  forces  acting 
between  the  atoms  in  a  molecule3.    First,  there  are  the 

forces  of  a+tra?tion  which  are  often  called  the  valence 
forces  and  which  fall  off  rapidly  as  the  distance  between 
the  atoms  becomes  larger.     Secondly,  there  are  forces  of 
repulsion  which  are  negligible  at  large  distances  but  which 

3.     Slater,  J.  C.  Introduction  to  Chemical  Physics  (New 
Trork$    McGraw-Hill,  1939)  p.  130 


increase  even  more  rapidly  than  the  attraction  at  small 
distances.     At  the  point  where  the  attraction  and  repul- 
sion forces  just  balance  is  the  equilibrium  position.  This 
is  the  normal  distance  of  separation  for  any  two  atoms, 
and  any  small  displacement  from  this  position  will  ?ive 
rise  to  forces  tending  to  restore  the  atoms  to  the  equili- 
brium position.     This  point  is  then  also  a  minimum  poten- 
tial energy  position,  and  the  amount  of  work  required  to 
pull  the  two  atoms  apart  is  the  energy  of  dissociation. 

When  two  atoms  approach  each  other,  there  is 
first  of  all  a  small  force  between  them  which  is  the  Van 
der  Waal's  attraction4".    This  force  arises  from  the  polar- 
ization of  each  atom  by  the  fluctuating  dipole  of  the  other 
and  the  overall  effect  is  rather  small.     As  the  atoms  ap- 
proach still  closer,  the  electronic  charge  of  one  is  at- 
tracted by  the  nuclear  charge  of  the  other.    This  is  es- 
sentially a  coulomb  attraction  even  thouph  each  atom  is 
electrically  neutral.     This  results  when  the  electronic 
charge  distribution  of  one  overlaps  the  electronic  charge 
distribution  of  the  other.     As  the  overlapping  increases 
and  the  nucleus  of  one  begins  to  penetrate  the  electronic 
charge  distribution  of  the  other,  strong  repulsion  sets 
up  between  the  two  nuclei. 

A  still  more  important  factor  must  be  considered, 
and  that  is  the  fact  that  when  atoms  be?in  to  overlap,  a 


4.     Ibid,  p.  368 
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distortion  of  charge  distribution  takes  place.    When  two 
atoms  overlap,  some  of  the  charge  shifts  from  the  region 
between  the  atoms  to  the  far  sides  of  the  nuclei.  This 
results  in  an  increase  of  energy  and  hence  a  repulsion  be- 
tween the  atoms.    However,  this  overlapping  and  distortion 
also  results  in  a  force  of  attraction  which  is  greater  than 
the  previous  forces  of  attraction.     This  is  the  exchange 
effect  which  refers  to  an  exchange  of  electrons  between  the 
two  atoms.    This  results  in  a  piling  up  of  electrons  between 
the  two  atoms  and  a  lowering  of  the  potential  energy,  and 
hence  a  binding  effect  is  produced.     The  electrons  involved 
in  this  interaction  are  the  outer  or  valence  electrons  and 
the  type  of  attraction  is  called  homopolar  valence  attraction. 
It  is  necessary,  however,  from  the  exclusion  orinciple  that 
the  electrons  entering  into  the  exchange  effect  have  oppo- 
site spins. 

When  more  than  two  atoms  are  involved,  the  number 
of  atoms  which  can  enter  into  the  attractive  exchange  type 
of  interaction  is  determined  by  the  number  of  valence  elec- 
trons.   Only  two  electrons  can  enter  into  an  electron-pair 
bond  and  if  there  is  only  one  valence  electron  on  the  atom 
involved,  then  there  will  be  repulsion  to  any  other  atom  that 
may  approach.     If  an  atom  has  a  number  of  valence  electrons, 
however,  then  it  can  form  a  bond  with  a  number  of  other  atoms 
depending  on  the  number  of  valence  electrons. 
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3.      The  Rotation  of  Molecules 

To  a  fairly  good  approximation,  the  rotation  and 
vibration  of  the  atoms  can  be  treated  separately  considering 
the  vibrations  as  if  the  molecules  were  not  rotating  and  the 
rotation  as  if  there  were  no  vibration^.    This  assumes  that 
the  amplitude  of  vibration  is  small  as  comoared  to  the  inter- 
nuclear  distances  and  that  the  forces  between  atoms  that  are 
induced  by  rotation  are  small  as  compared  to  the  forces 
between  atoms  which  give  rise  to  the  restorine  forces  re- 
sulting in  vibrations. 

The  far  infra-red  bands  correspond  to  a  change  of 
rotational  quantum  number  alone  and  they  will  have  an  in- 
tensity only  if  the  molecule  possesses  a  permanent  dipole 
moment  since  otherwise  the  average  value  of  the  dipole 
moment  will  be  zero  and  no  dipole  radiation  will  occur. 
Any  molecules  which  possess  electrical  symmetry  will  have 
no  permanent  dipole  moment,  and  hence  the  only  rotational 
spectra  which  can  be  observed  in  these  cases  are  those  that 
arise  due  to  the  formation  of  a  dipole  moment  because  of  un- 
symmetrical  vibration  of  the  atoms.     These  are  generally 
observed  as  fine  structure  on  the  vibrational  band  spectra. 
We  are  here  concerned  primarily  with  the  vibrational  spectra 
of  molecules. 

5.      Dennison,  David  M.,  Rev.  Mod.  Phy.  3,  280  (19?1) 


4.      Vibration  of  Molecules 

Y/hen  the  atoms  in  a  molecule  are  in  their  equi- 
librium position,  there  are  forces  actine  on  them,  as  previ- 
ously discussed,  which  will  tend  to  return  them  to  this 
position  if  they  are  displaced.    This  point  is  a  potential 
energy  minimum  and  any  displacement  tends  to  increase  the 
potential  energy.     If  the  force  function  tending  to  restore 
the  atoms  to  the  equilibrium  position  is  linear  in  this 
neighborhood,  then  for  small  displacements  the  atoms  will 
perform  simple  harmonic  oscillations  about  this  position. 
Dennison  has  shown  that  under  these  conditions  the  fre- 
quency of  radiation  absorbed  or  emitted  is  very  nearly  that 
of  the  classical  frequency  of  mechanical  vibrations.  Since 
infra-red  absorption  bands  are  concerned  with  transitions 
from  the  lowest  energy  states  to  those  immediately  above  it, 
it  is  possible  to  consider  only  vibrations  of  the  lowest 
energy  states  where  the  force  function  is  nearly  linear  and 
hence  classical  methods  may  be  used  for  the  determination 
of  the  vibrational  frequencies.    However,  to  determine  in- 
tensities and  selection  rules  quantum  mechanics  must  be  used. 

In  general,  the  normal  modes  of  vibration  of  a 
molecule  are  determined  as  a  function  of  the  mutual  dis- 
placements of  the  particles  from  their  equilibrium  positions. 
While  normal  vibrations  may  be  determined  by  considering  the 
forces  acting  on  the  oarticle,  in  practice  usually  the  con- 
servation of  energy  principle  is  utilized.    This  requires 
setting  up  the  kinetic  and  potential  energy  in  terms  of  the 
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displacement  coordinates.     The  type  of  potential  function 
that  is  used  is  determined  to  a  certain  extent  by  the  type 
of  molecule  which  is  being  considered. 

5>      Potential  Functions 

Usually  the  potential  function  is  developed  as  a 
power  series  in  the  displacement  coordinates  in  which  case 
all  linear  terms  are  absent  since  they  refer  to  the  rest 
position:  the  quadratic  terms  give  the  normal  vibrations, 
and  higher  order  terms  give  an  harmonic  corrections.  Several 
different  choices  of  the  coordinates  are  possible  and  the 
potential  function  will,  to  a  certain  extent,  depend  on  the 
choice  of  these  internal  coordinates. 

If  the  center  of  mass  and  angular  momentum  of  the 
molecule  are  not  considered  fixed,  there  will  be  y/  coordi- 
nates where  N  is  the  number  of  atoms  in  the  molecule.  How- 
ever, since  there  is  no  restoring  force  for  simple  trans- 
lation or  for  rotation  about  any  axis,  these  motions  give 
non-genuine  normal  vibrations.     Thus,  if  we  fix  the  center 
of  mass  and  set  the  angular  momentum  equal  to  zero,  the 
number  of  coordinates  will  be  reduced  by  five  if  it  is  a 
linear  molecule  or  by  six  if  it  is  non-linear.    This  will 
leave  either  3N-5  or  internal  coordinates.     One  choice 

of  coordinates  will  then  be  to  use  3^-6  internuclear  dis- 
tances and  deviations  from  these  3N-6  internuclear  distances 
as  internal  coordinates.    Such  coordinates  are  sometimes 
called  central  force  coordinates  since  they  work  particularly 
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well  with  a  central  force  system^.     In  this  type  of  force 
field,  the  assumption  is  made  that  the  force  acting  on  any 
one  atom  is  the  resultant  of  the  attractions  and  repulsions 
of  all  the  other  atoms  in  the  molecule  and  depends  only  on 
the  distance  between  atoms  and  acts  along  the  line  connec- 
ting the  atoms.    This  type  of  force  field  has  been  found  to 
work  fairly  well  for  molecules  composed  of  heavy  atoms. 

The  type  of  internal  coordinates  that  seem  to 
work  best  in  actual  practice  are  symmetry  coordinates 7,8,9. 
These  are  based  on  the  symmetry  properties  of  the  molecules. 
For  a  given  value  of  one  of  these  coordinates  there  corre- 
spond displacements  of  the  nuclei  in  agreement  with  one  of 
the  symmetry  types.     These  symmetry  types  will  be  discussed 
late"^  in  the  discussion  on  symmetry  and  group  theory. 

In  the  establishment  of  the  potential  energy  ex- 
pression, there  are  usually  more  potential  constants  than 
normal  frequencies.     Hence  only  in  a  few  cases  is  it  possible 
to  obtain  the  values  of  the  potential  constants  on  a  purely 
theoretical  basis.     One  way  to  solve  this  difficulty  is  to 
consider  an  isotope  effect^  where  one  or  several  atoms  are 

6.  Shaffer,  W.  H.  and  Newton,  R.R.,  J.  Chem.  Phy..  10,  405 
(1942) 

7.  Howard,  J.  B.  and  Wilson,  E.  B.  Jr.,  J.  Chem.  ?hy.  2, 
620  (1934) 

8.  Redlich,  0.  and  Tompa,  H. ,  J.  Chem.  Phy.  5,  529  (1937) 

9.  Rosenthal,  J.  E.  and  Murphy,  G.  M. ,  Rev.  Mod.  Phy.  8, 
317  (1936) 

10.  Herzberg,  Gerhardt,  Molecular  Spectra  and  Molecular 
Structure. .  II.  Infra-red  and  Raman  Spectra  of  Polyatomic 
Molecules  (rew  York:  D.  Van  Nostrand  Co.,  Inc.,  1945) 

p.  227  ff . 
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replaced  by  isotoplc  atoni3  for  which  the  force  constants  are 
the  same  but  the  frequencies  are  different.     This  gives 
additional  equations  in  the  same  constants  and  permits 
solving  for  them.    Another  solution  is  to  make  special  as- 
sumptions about  the  forces  acting  in  the  molecules  thus  re- 
ducing the  number  of  constants.     One  such  assumption  is  the 
central  force  field  assumption  which  has  already  been  men- 
tioned. 

Another  type  of  force  field  assumption  is  the 
valence  force-^.     This  assumes  that  there  is  a  strong  re- 
storing force  along  the  line  of  ever^  valence  bond  and  an- 
other opposing  the  change  of  the  angle  between  two  valence 
bonds.    This  type  of  force  field  appears  to  be  a  better  ap- 
proximation in  most  cases  than  the  central  force  field. 

Another  type  of  more  general  force  field  could 
be  mentioned  and  that  is  the  Urey-Bradley  field12.  This 
type  of  field  seems  to  apply  fairly  well  to  molecules  con- 
taining heavy  atoms  in  a  tetrahedral  form.    A  repulsive 
interaction  is  assumed  between  the  corner  atoms  of  the  mole- 
cule thus  adding  terms  of  the  order  l/rjn  to  the  potential 
energy  expression. 

6.      Mathematical  Methods 

(a)  Classical  -  Since,  as  has  been  indicated, 
the  vibrations  of  the  atoms  in  the  molecule  for  small  dis- 
placements can  b e  considered  on  the  basis  of  classical 

11.  Wilson,  E.  B..  Jr.,  J.  Chem.  Phy.  9,  76  (1941) 

12.  Urev,  H.  C.  and  Bradley,  C.  A.,  Phy.  Rev.  38,  1969 
(1931)  " 
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vibrations,  it  is  possible  to  use  classical  methods  to  deter- 
mine the  normal  frequencies  of  vibration.     The  mathematical 
methods  outlined  here  follow  in  major  part  the  development 
as  given  by  Herzberg^.     In  practice,  the  energy  principle 
is  generally  utilized  in  developing  normal  vibrations  so  that 
is  the  method  that  will  be  outlined  here. 

Consider  a  molecule  which  has  N  atoms.     If  a  part- 
icle is  displaced  from  rest,  there  will  be  a  restoring  force 
which  will  depend  only  on  the  components  of  the  displace- 
ment from  the  rest  position  x^,  y^,  z^.    The  particle  oscil- 
lating in  simple  harmonic  motion  about  its  rest  position  will 
have  a  kinetic  energy 

T±  =  1/2       v2  =  1/2        b±2  (1) 

and  potential  energy 

Vi  =  2n2^i2  mi  si2  =  V2  *±  s12  (2) 

where 

ki  =  4n2^i2  mi-  (3) 
If  we  let  x^,  y^,  z^,  x2,  3^2*  z2*         ^e  re-Dresentea* 
by  q-p  <\2"  ^31  ^41  ^  we  can  wrlte 

V  -  1/2    J=  kij  q,  qj  =  1/2  kn  qx2  ♦  l/2k22  q,2 

+  k12  ql  q2  +,,#  ^ 
where  the  force  constants  k^  ~  kji 

The  kinetic  energy  is  then  given  by 

T  mi  U2  +  y±2  +  z2)  (5) 

13.      Kerzberg,  Op.  Cit..  p.  72  ff. 


or  T  =  1/2  -f-}  «13  ,±  Qj  (6) 
where  =  o  for  i^j,  and  an=a22=a33=ml* 

a44*a55"a66"*m2#"  (7) 
We  can  now  transform  to  new  coordinates  ^2«.. 

a  linear  transformation 

^iccll  VC12  VC13  V'" 
q2=c2  ^l+c22  ^2+c23  ^3+"# 
q3=c31  ^l+c32  r?2+c33  ^3+*" 


(8) 


qi=cil  Vci2  r?2+ci3  ^3+--' 

By  appropriate  choice  of  the  coefficients  c^  both  V  and 
T  can  be  brought  to  a  simpler  form  in  terms  of  the  new  co- 
ordinates . 

V  -  1/2  (X^  r^2  +  X.2»?22+...+X.i^i2+..  A3nr?3n2)  (9) 

t  =  1/2  (n12^22+...+\2+...+n3n2)  (10) 

The       are  roots  of  the  secular  equation  ^ 


^"11  all^*  ^12  a12^*    ^"13  a13^v  *  *  * 

mm  mm  _ 

^21  a21^*  ^22  ^22^*    ^"23  ^23^*  •  •  • 


-  o  (11) 


The  total  energy  is  now 

H=T+V  -  l/2(»712+X.1^12)+l/2(ri22+X.2i722)  +  .<.  (12) 

It  is  thus  the  sum  of  3N  mutually  independent  terms  each 
of  which  has  the  form  of  a  simple  harmonic  oscillator.  The 


14.      Whittaker,  E.  T. ,  Analytical  Dynamics.  4th  Ed. 
(New  York,  Dover  Publications,  1944)  p.  181 


-  r 
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motion  of  the  entire  system  is  then  equivalent  to  the  super- 
position of  3N  simple  harmonic  motions  in  the  normal  co- 


ordinates. Thus! 

.0 


n±stl{  cos(2TTVit+(61) 


where  now 


2  2 


We  have  that 


k    =kn    k    =kn  k  =kn 

Kll  Kxx«  K22  Kyy          K12  Kxy' 


(13) 


(14) 


(15) 


Substituting  these  values  and  those  from  Eq,  (7)  in  Eq.  (11) 
gives 


XX  1 

k11 

Kxy 

k12 

*xx  ••• 

,1N 
Kxz 

k11 

kyy-mlX 

* 

k12 

Kyx  ••• 

k1N 

k11 

Kzx 

kn 
zy 

kzz"ml^  kzx  • 

k1N 

•*  *zz 

k21 

Kxx 

k21  k21 
Kxy  *xz 

kxx~m2^ 

k2N 

xz 

=  0 


(16) 


The  form  of  any  normal  vibration  is  obtained  by  putting  all 
other  t^'s  in  Eq.  (8)  equal  to  zero.    Then  for  ^ 


U)= 


(17) 


As  previously  discussed,  either  five  or  six  of  these  roots 
are  non-genuine  and  correspond  to  translation  and  rotation^ 

hence  there  are  actually  only  3^-5  or  3N-6  normal  vibrations. 
This  still  gives  a  large  secular  equation  to  solve  even  if  N 
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is  only  moderately  large.    A  considerable  amount  of  simpli- 
fication can  be  realized  if  the  molecule  possesses  any 
symmetry.     The  application  of  group  theory  to  polyatomic 

15 

molecules  is  discussed  in  detail  by  Rosenthal  and  Murphy  J 
and  others16'17. 

It  is  shown  in  the  references  cited  that  if  a 
symmetry  operation  is  carried  out  which  transforms  the  mole- 
cule into  a  form  indistinguishable  from  its  original  form, 
then  the  force  field  and  potential  energy  will  also  be  the 
same  as  before  and,  of  course,  the  normal  vibrations  will 
be  the  same.    However,  this  need  not  necessarily  be  the  case 
for  the  displacements  in  a  vibrating  molecule.    The  N  normal 
frequencies  may  then  be  divided  into  a  number  of  symmetry 
classes  each  of  which  will  contain  a  given  number  of  fre- 
quencies.   These  classes  or  irreducible  representations, 
as  they  are  called  in  group  theory,  will  all  be  along  the 
princioal  diagonal  of  the  secular  equation  and  can  be  treat- 
ed as  separate  equations.    These  symmetry  operations  aid  in 
assigning  the  normal  modes  of  vibration  to  certain  displace- 
ments of  the  atoms. 

The  actual  determination  of  the  frequencies  of 
vibration  now  requires  the  complete  solution  of  the  secular 
equation  which  can  usually  only  be  done  for  very  simple  cases 
if  the  force  constants  are  not  known.     In  order  to  reduce  the 

15.  Rosenthal,  J.E.  and  Murphy    Rev.  Mod.  Phy.t  8    317  (19?6) 

16.  Meister,  A.G.,  Cleveland,  F.F.  and  Murray,  Iff. J.,  Am.  J. 
Phy..  11,  239  (1943) 

17.  Herzberg,  Op.  Cit..  P.  1  ff.  and  p.  82  ff. 
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number  of  force  constants,  it  is  necessary  to  introduce  a 

force  field  assumption.    The  various  types  of  assumptions 

that  have  been  made  and  the  manner  in  which  the  calculations 

18 

are  then  carried  out  are  given  in  detail  in  Herzberg  • 

When  the  molecule  is  rather  large,  the  calculations 

become  complex  and  it  is  almost  impossible  to  carry  out 

algebraic  calculations  and  obtain  expressions  relating  the 

force  constants  to  the  frequencies.    The  computations  are 

then  best  carried  out  numerically.    An  important  factor  that 

is  of  great  aid  in  these  cases  is  the  fact  that  both  bond 

stretching  and  bending  force  constants  are  invariant  from 

one  molecule  to  another  if  the  bonds  are  in  similar  surround- 
19  20 

ings     "     •    Also  for  a  general  check  on  results  character- 
istic groups  of  atoms  have  vibrational  frequencies  that  are 
essentially  the  same  in  different  molecules  if  their  fre- 
quencies differ  sufficiently  from  the  vibrational  frequencies 

21 

of  the  rest  of  the  molecule     •    The  calculated  values  of  the 
frequencies  can  then  be  used  as  an  aid  in  assigning  the  ob- 
served bands  to  various  atomic  vibrations. 

(b)  Quantum  Mechanical  -  While  the  calculation  of 
the  normal  frequencies  can  generally  be  carried  out  entirely 
by  classical  methods,  the  frequencies  can  be  determined  by 


18.  Herzberg,  Op.  Cit..  Dp,  159-191. 

19.  Crawford,  B.L.,  Jr.,  J.  Chem.  Phy.  7,  555  (1939) 

?0.      Crawford,  B.L.,  Jr.,  ana  brinKiey,  SIR.,  Jr.,  J.  Chem. 

Phy.9T  69  (1941) 
21.      q^s5j'  E.B.Jr.,  and  Wells,  A.J.,  J.  Chem.  Phy.  9,  314 


quantum  mechanical  methods  and  the  selection  rules  are 
determined  by  quantum  treatment.     In  a  wave  mechanical  treat- 
ment, the  vibrational  motion  of  the  molecule  as  in  the 
classical  method  may  be  considered  as  a  superposition  of  3~ 

simple  harmonic  motions  to  a  good  approximation.  The 

22 

Schrodinger  equation  in  normal  coordinates  is 


where 


E  =  E1+E2+...+E3n.  (2) 


The  eigenf unctions  of  Eq.  (1)  are  the  eigenf unctions  of  the 
ordinary  simple  harmonic  oscillator 

where  N  .  is  a  normalization  constant,^  •*  ^ZH^'and 

VI  f\ 

(/5c7"^wc)  is  a  Hermite  polynomial  of  the  v^th  degree. 
The  total  vibrational  eigenfunction  is  then  the  product  of 
the  3^-5  or  harmonic  oscillator  eigenfu net ions . 

The  eigenvalues  of  the  ith  harmonic  oscillator 
are  given  by 

Ei=hvL  <vi+1/2)'  Tis0,l,2,.M  (4) 

where 

Vi=l/2n/5C^  (5) 

is  the  classical  frequency  of  vibration  of  the  ith  normal 


22.      See  Herzberg,  Op.  Cit.;  Pauling,  L.,  and  Wilson,  E.B. 
Jr.,  Introduction  to  Quantum  Mechanics  (Few  York: 
McGraw-Hill,  1935)  or  Eyring,  Henry,  Walter,  John,  and 
Kimball,  George  E.,  Quantum  Chemistry  (Few  York:  John 
Wiley  and  Sons,  1944) 
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vibration.    The  total  vibrational  energy  of  the  system  can 
then  have  the  values 

■(▼li  v2*  v3,,,)  =  h^1(v1+l/2)+h^2(v2+l/2)+...  (6) 

The  eigenfu.nct ions  are  either  even  or  odd  func- 
tions of  the  normal  coordinates  depending  on  whether  the 
energy  quantum  number  is  even  or  odd.     Thus  if  a  normal  vibra- 
tion is  antisymmetric  with  respect  to  a  symmetry  element, 
Vj^^)  changes  sign  for  odd  v^  but  remains  unchanged  for 
even  v^  if  the  particular  symmetry  operation  is  carried  out. 
Since  the  total  vibrational  eigenfunction  is  a  product  of 
the  harmonic  oscillator  eigenfunctions,  for  only  non- 
degenerate  vibrations  the  total  eigenfunction  will  be  sym- 
metric with  respect  to  a  given  symmetry  operation  if  the 
total  number  of  components  of  the^/AC^)  that  are  anti- 
symmetric is  even;  and  of  course,  it  will  be  odd  if  the 
total  number  of  antisymmetric  components  is  odd.    When  there 
are  degenerate  vibrations,  each  of  the  degenerate  eigen- 
functions goes  into  a  linear  combination  of  the  eigenfunctions 
belonging  to  the  same  energy  level. 

There  are  only  certain  combinations  of  symmetry 
properties  of  the  normal  vibrations  and  vibrational  eigen- 
functions possible  because  only  certain  combinations  of  sym- 
metry elements  occur  in  each  point  group  and  some  of  these 
are  necessary  consequences  of  others.    These  combinations 
are  the  symmetry  types  or  soecies  and  are  the  irreducible 
representations  of  group  theory.    For  example,  if  a  point 
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group  has  two  necessary  elements  of  symmetry  there  are  four 
possible  symmetry  types  ++,  — ,  +-,  -+.    These  species 
along  with  the  character  tables  have  been  worked  out  for  all 
of  the  different  point  groups^. 

Now,  since  the  vibrations  of  the  atoms  have  been 
considered  as  a  superposition  of  simple  harmonic  oscillators, 
the  selection  rule  to  a  good  approximation  is  given  by 

LV±  •  +  1  (7) 
for  each  normal  vibration  v^.    However,  the  presence  of  a 
fundamental  in  either  the  infra-red  or  Raman  spectra  depends 
on  the  change  of  the  dipole  moment  or  polarizability. 

The  dipole  moment  of  a  molecule  is  represented  by 

O  A 

the  matrix 

}  (8) 
where  M  is  a  vector  with  components 

!V=£<=ix>>  M^=Z4rtyi,M*sXei%  (9) 
where  e^  is  the  charge  on  the  ith  particle  with  coordinates 
xi*  ^i*  zi*  anc^  ^n  andVm  are  the  time  dependent  eigenf unctions 
in  the  states  m  and  n,  then 

%  *  V„  e^*W*  ,    %.MWe      T    .  (10) 
The  diagonal  elements  represent  the  permanent  dipole  moments 
in  the  states  n;  the  off-diagonal  elements  correspond  to  a 
transition  from  the  state  n  to  m.    The  transition  proba- 
bility is  proportional  to  the  square  of  the  time  dependent 

component  of  Eq.  (8).     The  vibrational  selection  rule  is 

23.  See  Herzberg,  Op.  Cit.    Rosenthal      Murphy,  Op.  Cit.. 
Kohlrausch,  K.  W.  F. ,  Per  Smekal-Raman-Sf f ekt  (Berlin; 
J.  Springer,  1931),  etc. 

24.  Herzberg,  Op.  Cit..  p.  251 


then  given  by  the  condition  that  at  least  one  of  the 
quantities 

^'^'^^  V%r'^M^Tfv*vC-/1*  (11) 

is  totally  symmetrical.     The  rule  may  also  be  stated  as:  a 
vibrational  transition  /wWVis  allowed  only  when  there  is 
at  least  one  component  of  the  dipole  moment  M  that  has  the 
same  species  as  the  product  tf^-'  t/V'« 

For  the  Raman  effect  the  intensity  depends  on  the 
induced  dipole  P  which  is  represented  by 

\        P*-4v  Jr  (12) 

where 

PjC  =  oiy>r  E  *-»-  ot^E^-r  o(XJ.Ej-  (13) 

Pfr*«tyl   Ey+^E^-KtyvEi.  (llf) 

P*  =  <*x  Ex  +rtz^  E^  ^-^aL«.E*  (15) 

where  the      are  the  components  of  the  polarizabilitjr  and  Hf 
is  the  electric  vector.    The  time  independent  part  of  Eq. 
(12)  is 

rpOjnm  =  ^  pO*^m  ^  (l6) 

o 

where  P    is  the  amplitude  of  P.    The  intensity  cf  a  Raman 
transition  is  proportional  to  the  square  of  JP0]11111. 

If  the  values  for  P  are  substituted  in  this 
equation,  we  get  the  components 

CP/r"  E\^w  V*y«Vf  *E)S*7  y>h  ytjtrtBlfaiWSJz  (17) 

P**£,fcVW  y^yfi  d-k  ^j^y,  v^Mf^«y*K^  (19) 
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Here  E°x,  E°y,  E°z  are  the  components  of  the  amplitude  of  the 

light  wave  and  the  integrals  are  the  matrix  elements  of  the 

six  components  of  the  polarizability  tensor.    The  diagonal 

elements  correspond  to  Rayleigh  scattering  while  the  off 

diagonal  elements  correspond  to  Raman  scattering.     The  Raman 

selection  rule  is  then  given  by  the  fact  that  at  least  one 

of  the  six  integrals  of  polarizability  must  be  different  for 

a  Raman  transition  to  occur  or  again  a  Raman  transition  is 

allowed  if  at  least  one  of  the  six  products 

cLtx  ^v'VC")  dry  Vv-'V^, ... 
is  totally  symmetrical. 

C.     THE  NORMAL  VIBRATIONS  OF  1,  3 -Butadiene 

1.    The  Configuration  of  lt  ^ .  Butadiene 
The  normal  configuration  for  butadiene  could  con- 
ceivably be  either  in  the  trans-  or  cis-  form  and  no  final 
conclusion  has  yet  been  reached  as  to  which  form  is  the 
lowest  energy  state.    Mulliken  has  shown  '  that  it  is  possible 
to  have  isomerism  about  a  single  bond  as  might  be  the  case 
with  butadiene  and  his  calculations  with  molecular  orbitals 

seem  to  indicate  that  the  trans-form  is  the  lowest  energy 

26 

form.    The  work  of  Carr,  Pickett  and  Stucklen      in  the  ultra- 
violet absorption  snectrum  of  butadiene  seems  to  indicate 

that  the  trans-form  predominates  at  room  temperature  at  least. 

27 

Rasmussen,  Tunnicliff  and  Brattain      have  obtained 

25.  Mulliken,  R.S.,  Rev.  Mod.  Phy.  14,  269  (194-2) 

26.  Carr,  E.P.     Pickett,  L.W..  and  Stucklen,  H.    Rev.  Mod. 
Phy.  14,  260  (1942) 

27.  Rasmussen,  R.S.,  Tunnicliff,  D.D.,  and  Brattain,  R.R.. 
J.  Chem.  Phy.  11,  432  (194^) 
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the  infra-red  spectrum  of  butadiene  and  compared  it  with  the 

Raman  spectrum  taken  with  the  butadiene  in  the  liquid  state 

28 

at  the  temperature  of  dry  ice  by  Bradacs  and  Kahovec  . 
They  found  a  coincidence  between  six  frequencies  in  the  Raman 
and  infra-red.    This  seems  to  indicate  that,  since  no  coinci- 
dences should  occur  if  the  butadiene  was  in  the  trans-form 
due  to  its  center  of  symmetry,  the  cis-form  is  predominant  at 
the  temperature  of  dry  ice.     Thus  it  would  appear  that  the  cis- 
form  is  the  lower  energy  form  but  at  room  temperature  the  trans- 
form predominates. 

The  cis-form  (Fig.  1)  has  been  selected  as  the  con- 

29 

figuration  for  the  present  work.    The  work  of  Mulliken  '  and 

28 

Bradac  and  Kahovec      indicates  that  all  the  atoms  are  in  a 
plane  and  that  all  bond  angles  both  H-C-H  and  C-C-H  are  very 
nearly  120° .    The  carbon-hydrogen  bond  distances  are  fairly 
well  known  from  electron  diffraction  studies-^    and  were  taken 
to  be  1.09  angstroms.    The  carbon-carbon  single  bond  and  car- 
bon-carbon double  bond  distances  were  taken  from  the  work  of 
Mulliken  who  based  his  conclusions  on  the  calculations  of 
Lennard-Jones-'    and  PenneyJ  •    The  C-C  distance    was  chosen 
as  1#44A.  and  the  C=C  distance  as  1.35A. 

28.  Bradacs,  K. ,  and  Kahovec,  L.,  Zeits.  fr.  Physik.  Chemie 
B48,  63  (1941) 

29.  Mulliken,  R.  S.,  J.  Chem.  Phy. ,  7,  121  (1939) 

30.  Pauling,  L.,  Nature  of  the  Chemical  Bond  (Cornell  Univ- 
ersity Press,  1939) 

31.  Lennard-Jones,  J.  E. ,  Proc.  Roy.  Soc.  158A,  293  (19^7) 

32.  Penney,  W.  G. ,  Proc.  Roy.  Soc.l5bTT"3l8  (1937) 


Fig.  1 

The  following  notation  has  been  adopted  in  the  work  that  fol- 
lows.   The  carbon  atoms  have  been  labelled  with  numbers  from 
1  to  4  while  the  hydrogen  atoms  have  been  lettered  with  Greek 
letters  as  shown  in  Fig.  1.    The  masses  of  the  carbon  atoms 
have  been  denoted  by  M  and  the  hydrogen  atoms  by  ra.  The 
C-C-C  angles  have  been  denoted  by  6  and  the  C-C-H  and  H-C-H 
angles  by         The  principle  axis  of  symmetry  is  chosen  as  the 
Z  axis. 

With  the  configuration  that  was  decided  upon  the 

point  group  for  the  molecule  is  C?    with  the  character  table 

v 

as  given  in  Table  1J  . 


33.    See  Rosenthal  and  Murphy,  Rev.  Mod.  Phys..  8,  317  (1936) 


Table  1 


£ 

9, 

/ 

1 

1 

( 

3* 

/ 

-1 

-1 

1 

ft* 

/ 

/ 

-/ 

-/ 

.  ..  &i 

I 

/ 

/ 

f 

.  <Pr 

0 

77" 

0 

6 

10 

0 

10 

0 

3 

-1 

/ 

1 

30 

0 

/o 

0 

2¥ 

3l 

JO 

0 

where  J0  is  the  angle  of  rotation.  Up  is  the  number  of  particl 
left  unchanged  by  any  operation,  <XR=UR(+l+2cos0^)  with  +  for 
proper  rotation  and  -  for  improper  rotation  and  5=(U^-2)x 
(l+2cosj8)  or  6  =Ug  (-l+2cosQ) .    From  the  selection  rules  the 
frequencies  in  species  A^,        and  B2  are  infra-red  active; 
A^,  A2,        and  B2  are  Raman  active.    The  Raman  lines  in 
A^  are  polarized;  those  in  A2,  B^  and  B2  are  depolarized. 
From  the  relation 

a(k)  =  \  (^(Uc-2)(l+2cos0)X(k)+  f  Us(-l+2cos0)X(k)}  (1) 

the  number  of  frequencies  in  each  species  may  be  calculated^4- 
Applying  this  gives  the  following  results 

9  (2) 

3  (3) 

4  (4) 
8  (5) 

When  the  symmetry  operations  are  applied  to  the 
molecule  together  with  the  character  sign  (either  symmetric 

34.    Rosenthal  and  Murphy,  Loc.  Cit. 


avx;  =  ±(24+2+10)  ■ 

a(2)  =  1(24-2-10)  = 

a(3)  =  1(24+2-10)  = 

a(4)  =  1(24-2+10)  = 


or  ant i symmetric)  from  the  character  table  and  the  conditions 
that  the  center  of  mass  shall  be  fixed  and  the  angular  momen- 
tum shall  be  zero,  the  following  cartesian  coordinates  are 
the  only  ones  that  remain  for  each  species: 

b2:  y,>  y-j^y-*  . 

The  only  vibrations  which  will  be  considered  are  in-plane 
vibrations;  hence,  species  A2  and  B2  will  give  no  frequen- 
cies since  these  vibrations  are  out  of  the  plane. 
2.    The  Potential  Function 

The  force  field  which  was  chosen  for  this  deter- 
mination was  the  valence  type  force  field.    Thus  the  changes 
in  bond  length  and  changes  in  the  bond  angles  were  chosen 
as  the  coordinates  for  the  potential  function.     The  kinetic 
energy  expression  can  be  written  directly  in  cartesian  co- 
ordinates as 

The  potential  energy  in  symmetry  coordinates  is 


(2) 
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The  values  of  |Ar12l2,  |Ar2^|2,...  are  determined  by  taking 
the  vector  difference  of  the  distances  between  two  atoms  in 
the  rest  position  and  when  they  have  been  given  a  small  arbi- 


trary  displacement.     These  are 

f J r y v - ^) r y v - y < •  - y ) J +  r* v - * r) A  -  K, 1 1 3 

- 


—  COS 
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\uTl\       UH^d  •  /    J  -l  i  ■   i—  


1^3 1 


rn3        *  -   —  —  


J. 


iLr.^':r'fr,''ii,itJ>'t^i-/^---ivi  -i 


aft 
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Calculation  of  the  Vibrational  Freouer.cies 
Now  from  the  application  of  group  theory  wp  have 
that  for  species  A-,  all  X's  are  zero  and 

«  )j  -  St** 

and  ^=  2=3-  •  rT5J  2r'  '  (3) 

For  species        ail  X's  are  zero  and 

where  the  constants  have  the  following  values: 

?  -    J^-jJ+C  f    -     la-'jJ+C  y  _  

~^7f     ^  5  ~ry>±) 


These  substitutions  are  made  in  the  expressions  given  for 
the  symmetry  coordinates  and  these  are  then  expanded  in  a 
Taylor's  series  leaving  only  second  order  terms.    The  final 
expressions  are  then  substituted  into  the  potential  energy 
expression  given  previously  and  the  coefficients  of  the 
cartesian  coordinates  are  grouped  together  giving  the  constants 
b^j  in  the  secular  equation 

I  bij  "  ^*1J  I  =  °*  (5) 
The  constants         are  similarly  obtained  by  applying  the 
symmetry  conditions  above  to  the  kinetic  energy  expression 

and  grouping  the  coefficients  of  the  cartesian  coordinates. 


The  constants  can  now  be  evaluated  numerically  and 
substituted  into  the  secular  equation.    The  values  for  the 
constants  that  were  used  throughout  this  work  are  given  in 
Table  II.     The  masses  of  the  atoms  were  taken  from  Birge'sjy 
table  of  values  of  physical  constants.    The  values  for  the 
force  constants  were  transferred  from  the  work  by  Wilson 
and  'Veils-     on  propylene  and  the  paper  by  Crawford  and  Brinkley 
on  force  constants  in  organic  molecules. 

Table  II 


This  gives  two  secular  equations,  one  of  eigth 
order  and  one  of  ninth  order.     The  expansion  of  secular 
equations  of  this  size  is  a  laborious  task,  and  usually  approxi 
mations  of  various  kinds  or  iteration  schemes  are  used  for  sol- 
ving them.    A  very  powerful  method  for  an  exact  expansion  of 
determinants  of  this  type  which  was  used  in  this  work  will  now 
be  described. 

35.  Birge,  R.T.,  Rev.  Mod.  Phys.  13,  233  (1941) 

36.  Wilson,  E.B.,  Jr.  and  '.Veils,  A.J..  J.  Chem.  Phy.  9,  314 
(1941) 

37.  Crawford,  B.L.,  Jr.  and  Brinklev.  S.R.,  Jr.,  J.  Chem. 
Phys.    9,  69,  (1941) 
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(a)    An  Exact  Method  for  Expanding  the  Secular 
Equation  -  The  method  here  described  is  due 
to  J.  G.  Bryan-5    who  devised  the  method  to  permit  an  exact 
expansion  of  the  characteristic  equation  and  latent  vectors 
of  a  matrix.    The  characteristic  equation  is  defined  by 

|a  -  m|  -  o  (6) 

where  I  is  the  unit  matrix.     The  process  of  expanding  this 
equation  consists  of  building  an  auxiliary  matrix  and  row 
vector  from  the  original  matrix  A. 

A  vector  f^  is  built  as  follows 

fl  =  "  1    all  *  (7) 
This  vector  is  then  used  as  the  basis  for  building  the 
auxiliary  matrix  l^*    A  matrix  A2  is  written  down  first  which 
is  composed  of  the  first  two  elements  of  the  first  two  rows 
of  A.     This  is  divided  into  two  parts  with  the  last  row 
being  indicated  by  a2  and  the  first  row  by      .  Thus 

h--rl-ll-   t  ^ 

a21    a22  a2 

The  auxiliary  matrix  B2  is  then  built  by  first  writing  down 
the  vector  f^.    The  first  row  is  built  by  multiplying  oi  in 
A2  times  a  null  column  for  the  first  term,  and  then  multi- 
plying o£  times  this  first  column  so  built  for  the  second 
term.  Thus 

B2  =  0       -a12  (9) 
fl  "  -1  all 

38.    Bryan,  J.  G. ,  A  Method  for  the  Exact  Determination  of 
the  Characteristic  Equation  and  Latent  Vectors  of  a 
Matrix  with  Applications  to  the  Discriminant  Function 
for  More  than  Two  Groups.    Harvard  Graduate  School  of 
Education  Dissertation  (1950) 
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The  process  is  then  repeated  by  writing  down  the 
first  three  terms  from  the  first  three  rows  of  the  original 
matrix  for  A^.    That  is 

all      a12  a13 
A^    =  a21      a22      a23  (10) 

a31     a32     a33  a3 

where  the  last  row  is  now  designated  by  a^  and  the  remaining 
rows  are  designated  by<*.    The  vector  f2  is  now  built  by 
multiplying  a2  by  the  null  column  in  B2  and  f^  and  subtract- 
ing the  next  number  in  f^,  that  is,  -1  for  the  first  term  of 
fp»    The  second  term  is  then  obtained  by  multiplying  a2  by 
the  first  column  of  B2  and  f^  and  subtracting  the  next  term 
in  f^.    The  last  term  is  obtained  by  multiplying  a2  times  the 
last  column  in  B2  and  f^.  Thus 

f2  ~  1  —  ^^22    all^     all  a22  (^-^-) 
The  matrix        is  then  built  in  the  same  manner  as  B2  was. 
First  f2  is  written  down  and  then  the  first  term  of  the  first 
column  of  B^  is  obtained  by  multiplying  the  first  row  of  °£  in 
A^  by  the  null  column  in  B^  and  f2»    The  second  term  in  the 
column  is  similarly  obtained  by  multiplying  the  second  row  in 
by  the  null  column.    Thus  there  will  always  be  a  column  of 
zeros  as  the  first  column  in  the  auxiliary  matrix.     The  first 
and  second  terms  of  the  second  column  in  B^  are  then  obtained 
similarly  by  multiplying  this  first  column  of    B^  and  f2  by 

the  first  row  of  oC  in  A^  for  the  first  term  and  by  the  second 
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row  of  oL  for  the  second  term.  The  last  column  in  is  simi- 
larly obtained.  The  vector  is  obtained  in  the  same  manner 
as  ?2         the  entire  process  is  continued  until  the  entire 

matrix  has  been  written  down  as  A  .     The  coefficients  of  the 

n 

expanded  characteristic  equation  will  then  be  given  by  the 

terms  in  f    which  is  obtained  in  the  manner  described  above 
n 

by  multiplying  aR  by  the  columns  of  Bn  and  fR_^  and  subtract- 
ing the  next  term  in,fn_^. 

There  are  several  advantages  to  this  method  besides 
the  obvious  advantage  of  permitting  an  exact  expansion  of  the 
characteristic  equation  with  even  less  labor  than  is  involved 
in  previous  approximation  methods.     One  of  these  is  the  fact 
that  the  method  is  readily  adaptable  to  use  by  I.B.M.  com- 
puting machines  or  hand  commuting  machines.     The  other  is 
the  fact  that  there  is  a  continuous  check  on  the  work  that 
has  been  performed.     The  check  is  obtained  by  summing  the 
columns  of  oi  and  multiplying  this  row  of  sums  times  the  col- 
umns of  B£  and  f-j^  that  is  being  built  by  this  particular oC . 
The  product  of  this  row  of  sums  of      in  some        times  a  given 
column  in  B^  and  f must  equal  the  sum  of  the  terms  in  the 
next  column  of  B^.    Furthermore,  once  B^  has  been  completed 
the  product  of  any  row  in oi times  the  last  column  in  B^  and 
^i-1  must  eo-ual  zero.    An  example  will  probably  serve  best 
to  clarify  the  operations  involved.     Consider  the  matrix 


2-\      -1  -4 
-2  3  1-A. 


=  0. 


This  can  be  written  as 


2    -1  -4 

4  3  2  -X. 
-2      3  1 


10  0 
0  10 
0    0  1 


=  0 


where  the  first  matrix  is  the  matrix  A  in  our  previous  dis- 
cussion. Then 


and 


Sum 


Sum 
Check 


fl  = 


A2 


B2 
fl 


-1 


2  ^1 

4-  3 

2  -1 


■1  2 

0  1 
0  1 


an 


Sura 


Sum 
Check 


"3 

(«0 

»s 

f2 


f3  = 


2 
4 


-1 


-4 

2 


?  3 

1 

a 

6  2 

-2 

0  -4 

10 

0 

0  2 

-20 

0 

1  -5 

10 

0  -2 

-10 

0  -2 

-10 

aL 

3 


-1     6    -1  -70 


Thus  the  expanded  characteristic  equation  in 

-  6\2  +  \  +  70  =  0. 
Now  the  secular  equation  considered  in  this  work 
is  in  general  not  in  the  form  of  the  characteristic  equation 
but  by  proper  addition  and  subtraction  of  rows  and  columns 
and  then  division  of  each  column  by  the  coefficient  of  \ 
in  that  column  after  there  are  X.'s  only  on  the  principal 
diagonal  it  can  be  put  into  the  proper  form.     Then  the  pro- 


r 


-  < 


r 


cess  described  can  be  used  to  obtain  the  expanded  secular 
equation.    As  a  further  example,  the  entire    solution  of  th 
eighth  order  secular  equation  is  presented.    The  original 


secular  equat 

2585.223001 

1566.097331 
285.972608 

211.?39583 
-I26I.519224 

-749.75037I 
-63.368231 
-I88.78209O 


ion  is  as  follows 
1. 

-  \919.445462 

-  X53^. 133809 

-  X.IOO.388702 

-  x  75.544711 

+  X.437. 160272 
+  X.225. 631124 
+  X  22.^70756 
+  X  65.401606 


1566.097331 
1031.210323 

172.922857 
132.222504 
-708.579705 
-413 . 541948 
-35.941267 
-106.43601S 


2. 

-  X 533. 133809 

-  M89.  394129 

-  X  55.872969 

-  X  44. 774-624- 
+  X.195- 288892 
+  MOO.  794274 
+  X  9.997498 
+  X  29.21^303 


4. 


285.972608 

172.922857 
34.489040 
24.163362 

■141.764534 

-85.888794 
-9.308504 

-21.970341 


-I26I.519244 
-708.579705 
-141.764534 
-105.753050 
695.144505 
397.096281 
30.793178 
99.841587 


■63.358231 
-35.941267 
-9.308504 
-5.412186 
30.793178 
20.538526 
5.676052 
5.409625 


x.100. 388702 

X  55.872969 
X  17.123452 

x.  B. 431028 
\  53.115492 
x.  27.414450 
x  2.718074 
x.  7.946372 


5. 

+  A.437. 160272 
+  X19 5. 288892 
+  x  53.115492 
+  x  36.714398 

-  X.308. 527845 

-  x.148.953422 

-  X  14.768355 

-  X  43.175750 

7. 

+  x  22.370756 

+  X.  9.993498 

+  X  2.718073 
+  x  1.878782 

-  x  14.768355 

-  x  7.622377 

-  X  2.429588 

-  x  2.209429 


211.739583 
132.222504 
24.163362 

22.942814 
■10  5.7530  50 
-62.626083 

-5.412186 
-15.822713 


X  75.544711 
X  44.774624 

x  8.431028 
x  8.018383 
x  36.71439B 
X  18.949368 
x.  1.878782 
x  5.492679 


6. 


•749.750371  + 
■413.541948  + 
-85.888794  + 
-62.626087  + 
397.096281  - 
254.216094  - 
20.538526  - 
60.830535  - 


X.225. 631124 

x.100. 794274 

X  27.414450 

x  18.949368 

XI 48. 953  422 

X.  96.809801 
x  7.622377 

X  22.284260 


8. 

■188.782090  +  X 
■106.436015  +  x 

-21.970741  +  X 

-15.822713  +  X. 
99.841587  -  X 
60.830535  -  X. 
5.409625  -  X 
16.397630  -  X. 


65.401606 
29.216303 
7.946372 
5.492679 
43.175750 
22.284260 
2.209429 
8.133182 


r  r 


After  adding  and  subtracting  rows  and  columns 
until  the  terms  in  \  remained  only  on  the  diagonal  and  then 
dividing  each  column  by  the  coefficient  of  \  in  the  column, 


the  following  matrix  resulted: 


1. 

2. 

3. 

4. 

2.067366-X. 
.260648 

1.056297 
.89^082 

-.212079 
-.170916 
.475468 
-.165267 

.297095 
1.510432-X. 
.259149 
.850230 
.174022 
-.274173 
-.932967 
.206441 

.631730 

.135973 
2.419432-x. 

-.739404 

.121093 
.113293 

-.026247 
.0334-35 

.283279 
.236602 
-.392156 
2.^11303->- 
-.112740 
-1.208686 
-.595243 
.017732 

5. 

6. 

7. 

8. 

-.33760^ 

.243037 
.322316 

-.565803 
.829401->. 

-.94634O 
-I.309634 
.203969 

-.044726 
-.062945 
.049572 

-.997172 

-.155567 

1.900935-*. 

1.056800 

-.007546 

.641858 
-1.104950 
-.059245 
-2.533311 
-I.IIO603 
5.451680 
6.498104-K 
-.375795 

-2.341203 

2.565703 
.791970 
.791960 

1.815131 

-.408514 
-3.943537 

2.082605-^ 

This  can  now  be  expanded  by  the  method  outlined  above.  Then 
f,  ■  -1  2.067366 


k0    2.067366  .297095  cL 

*  .26064«             1.510432  a2 

B0  0  -.297095 

f^  -1  2.067366 

2.067366  .297095  .631740  , 

A,       .260648  1.510432  .13597^ 

5  1.056297              .259149          2.419432  a 


Sum  (^)2.328014  1.807527  .767703 
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B. 


0 

0_ 

1 


.631730 

■3.577798 


-.913788 
-.116447 


3.045179 


Sum  (B-)0 
Check    5  0 


.76770^ 
.767703 


■1.0^0235 
•1.030235 


Sum  (oO 


2.067366 
.260648 
1.056297 


".8930«2 
3.384311 


.297095 
1.510432 

*  259149 
.«50230 


.631740 

.135973 
2.419432 


2.066676       3.18713  5 


.283279 
.236602 
-.^92156 


-.739404  2.U1303 


.127725 


a. 


E 


4 


0 
0 
0_ 

-1 


-.283279 
-.236602 

,-3?21?6 
5.997230 


1.290691 
1.041072 
-I.76 


-10.991 


-1.25224- 

3327* 


-.9 


1.679522 
6.372195 


Sum(BA)  0 
Check^  0 


Sum  (ot) 


B, 


0 
0 
0 
0_ 

1 


Sum(B^)  0 
Check'  0 


.127725 
.127725 


.568166 
. 568166 


-.505997 
-.505996 


2.067366  .297095  .631730  .283279 

.260648  1.510432  .13597^  .236602 

1.C56297  .259149  2.419432  -.?92156 

.893082  .850230  -.739404  2.31130^ 

-.212079  .174022  .121093  -.112740 

4.277393  2.916906  2.447731  2.439028 


-.337603 
.2430^7 
•3223i6 

-.^65803 
.829401 

-.338053 


-.3^7603  2.222578    -4.468156  2.593151 

.243037  -I.830229      4.132162  -2.446678 

.^22316  -1.969895     3.680384  -2.076983 

-.565803  3.060059    -4.686908  2.066986 

-8. 308533  24.116182  -28.452099  11.574^ 


2.067366 
.260648 

1.056297 
.893082 

-.212079 
-.170916 


-.338053 
-.338053 

.297095 
1.510432 

.259149  ; 
.850230  • 
.174022 
-.274173 


1.482514  -1.342518  .n6476 
1.482513   -1.342518  .136477 


.631730 

.135973 
.419432 

.739404 
.121093 
.113293 


.283279 
.236602 
-.392156 
2.311303 
-.112740 
-1.208686 


-.33760^ 

.243037 
.322^16 

-.565&P 
.829401 

-.946340 


-.044726 
-.062945 
.049572 
-.997172 
-.155567 
1.900935 


Sum(o04. 065314    3.090928    2.568824    2.326288      .491348  -1.210838 


r 


i 
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0 
0 
0 
0 
0 


•044726 
.062945 
-•049572 
.997172 
.155567 


-.098896 
-.201454 
.055700 
-6.765270 
-1.409489 


is    -1    9.137934  -30.790576 


Sum®  0 
Check  0 


1.210838 
1.210878 

2.067366 
.260648 

1.056297 
.893082 

-.212079 
-.170916 
.475468 


-8.419769 
-8.419369 

.297095 
1.510472 

.2S9149 
.850270 
.174022 
-.274173 


-.932967 


-.292578 
-.727598 
.650466 

15.563708 
4.376779 

47.080691 

19.970^37 
19.970337 

.631730 

.13  5973 
2.419432 

-.7^9404 

.121093 
.117297 

-.026247 


.53  4401 
1.299894 
-1.179141 
•14.4720^4 

-5*365307. 

■32.531812 

■19.1821«7 
■19.182188 

.287279 
.276602 
-.792156 
2.  71170-? 
-.112740 
-1.208686 


-.087166 
-.738005 
.381805 
4.480598 
2.212552 
8.179528 

6.249784 
6.249783 

.777603 
.247037 
.^22716 
-.565807 
.829401 
-•246740 


-.595243  -1.709634 


.044726 
-.062945 

.049572 
-.997172 
-.155567 
1.90093  5 
1.0  56800 


.641858 
■1.104950 
-.0  59245 

■2.533311 
■1.110607 
5.451680 
6.498104 


Sum(^)  7.894398      2.816755      2.682117        1.117602  -.454992 


B, 


0  .641858 

0  -1.104950 

0  -.059245 
0  -2.533311 

0  -1.110607 
5.451680 


-6.710656 
9.475230 
1.808043 

16.979335 
10.440  585 
•45.517575 


27.429747 
-29.459247 
-8.579664 
-^9.051691 
-34.874754 
176.844687 


{,     1-11.038869    46.778189  -96.022237 


.690097  1.285429 

-32.789724 

40.224360 

12.915582 

38.467198 

49  .  730^19 
■180.0824-91 

99.289523 


18.217245 
-23.221188 

-6.828934 
-15.820024 
-29.828843 
101.452147 


Sum(B7) 1.28 5429  -13.525038  48.349074 
Check    1.285429  -13.525037  48.349074 


-47. 992100 
-71.534756 
43.966799 
-71.534755 

43.966799 


-3.312372 
4.395140 
1.160792 
2.298122 
5.870528 

■18.95428 
8.22774 


-8.582073 


-8.582073 


a.. 


k8 


2.067366 
.260648 

1.056297 
.893082 

-.212079 
-.170916 
.475468 
-.165267 


297095 
510432 

259149 
850230 
174022 
274173 
932967 
206441 


.631730 

.135973 
2.419432 
.739404 

.121093 
.113293 

-.026247 
.033435 

-.044726 
-.062945 
.049572 
-.997172 
-.155567 
1.900935 
1.056800 
-.007546 


.283279 
.236602 
-.^92156 
2.311303 
-.112740 
■1.208686 


.017732 

.641858 

■1.104950 

-.059245 

•2.5333H 
■1.110603 

5.451680 

6.498104 

-.375795 


Sum(^)     4.369866  1 


E 


8 


0 

2.341203 

-34.578691 

0 

2.56570^ 

36.6IO239 

0 

-.79197c 

I3.292887 

0 

-.791960 

3.182598 

0 

-1.815131 

24.933551 

0 

.408514 

17.999896 

0 

3.943537 

-36.724235 

1. 

17.536973 

-108.4485^3 

883788  2.655870  .522359 
1.746897  7.783533 


169.794569 
-183.896618 

-70.908184 

-27.136810 
-123.227256 
-151.753276 

122.560845 

316.030387 

-170.460830 
163.557080 
93.02Q443 
180.066682 
184.716806 

241.720995 
-30. 148225 
-133. 881886 


-.337603 
.2430^7 
.322316 

-.565703 
.829401 
.946340 
-I.309634 
.203969 

-2.341203 
2.565703 
.791970 
.791960 
1.815131 
-.408514 

V?43^7 

2.082605 

-I.764626 

-.728490 

•368.924785 
405.450890 
166.565841 
145.839614 
285.^84102 
413.366798 
•178.539182 


377.377522 
■394.220247 
•188.719756 
•265.155967 
.331.574786 

•47D.760675 
115.150376 

361. 391323 


Sum  (Bg)  0  .728490        24.716245  -264.566730 

662.481951 
■264.566730 
662.481952 
.509.361764 
624.402846 


Check       0  .728490 


fg        1  -19.619578 


■1157.903  533 
24.716246 
•1157.903534 
142.158734 
■IO67.H6658 


-470.918704 

26.887299 
-22.645782 
-15.931958 
-40.977042 
-38.351770 
-47.787798 
2.334004 

18.698157 

869.143278 

-136.473047 
869.143278 

-136.473049 
990.160760 

-182.009395 


■s 


+  20.223912 

These  are  the  coefficients  of  the  expanded  secular  equation. 
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(b)    The  Determination  of  the  Formal  Freouencies  - 
The  solution  of  the  expanded  equation,  which  was  accomplish- 
ed by  Horner's  method,  yields  the  following  eight  roots: 

. 283  ?85xl029 ,    . 53  500 5xl029 ,      . 691997xl029 ,    1 . 0625l4xl029 
2.091375xl029,  3.10l65lxl029,    3.244094xl029,  8.609156xl029 

2  2 

But  these  are  equal  to  4rr         where  the        are  the  frequen- 
cies of  vibration.    Then  solving  for  the  wave  numbers  gives 
£    =  894  cm"1    ,  1      =    1228  cm"1  , 

J    =1397  cm"1    ,  J     =    1731cm"1  , 

A3  »  N 

-    =2428  cm"1    ,  |      =    2957  cm"1  , 


5 

h  =3°24  Cm"X     '  L    -    4926  cm"1 


1  .-1  1 


•8 

The  expansion  of  the  ninth  order  secular  equation 

and  the  solution  for  the  wave  numbers  by  the  same  method  as 

described  above  gives  the  following  nine  values: 

I    =434  cm"1  ,       I    =  1182  cm"1  ,      J      =1324  cm"1 
**!  2 


*3 

J    =  1642  cm"1  ,      J    =  1642  cm"1  ,      i      =    2062  cm"1  , 
*4  K5  H 

Z?  =  2846  cm      ,      ^  =  2992  cm-l         1^    =    2994  cm"l  ( 

The  values  for  2^  and  ^9  indicate  degeneracy;  how- 

ever, these  are  not  first  order  degeneracies.    The  pair  of 
roots  for  1/^  and  2^  in  the  solution  of  the  equation  turned 
out  to  be  just  barely  imaginary;  that  is,  the  curve  just 
barely  missed  crossing  the  axis.    A  recheck  of  all  of  the 


work  showed  that  there  did  not  seem  to  be  an  error  which 
indicates  that  one  or  possibly  more  of  the  potential 
constants  should  probably  be  adjusted  slightly.  However, 
it  may  also  be  the  case  that  these  two  pair  of  roots  may  be 
accidental  degeneracies.    The  general  grouping  of  the  lines 
seems  to  check  fairly  well  with  observed  spectra,  as  will 
be  demonstrated  later,  so  that  these  may  very  well  be  acci- 
dental degeneracies  produced  by  interactions  between  vibra- 
tions. 

4.    Comparison  of  Results  with  Observed  Spectra 
The  observed  frequencies  for  butadiene  that  are 
listed  in  Table  III  were  obtained  from  results  published  by 
Rasmussen,  Turmicliff  and  Brattain^^,  Creitz  and  Smith40, 
Richards  and  TTielson    ,  and  The  American  Petroleum  Institute 
The  calculated  frequencies  are  listed  in  the  adjacent  column 
with  matching  between  observed  and  calculated  frequencies 
indicated.    A  comparison  of  observed  and  calculated  spectra 
is  also  made  in  Fig.  2.    As  is  apparent  from  the  table  and 
figure,  there  seems  to  be  fairly  good  matching  of  the  general 
distribution  of  the  calculated  lines  with  the  observed. 

So  long  as  interaction  between  the  end  groups  of 
atoms  in  the  chain  is  not  taken  into  consideration,  as  was 

Rasmussen,  Turmicliff,  and  Brattain,  Loc.  Cit. 

40.  Creitz,  E.  C.  and  Smith,  F.  A.,  J.  Res.  Nat.  Bu.  Stand. 
43,  365  (1949) 

41.  Richards,  Charles  M.,  and  Nielsen,  J.  Rtfrd,  J.  Op.  Soc. 
Am.  40    433  (1950) 

42.  Am.  Petroleum  Inst.  Res.  Proj.  44  at  N.B.S.  Cat.  of  Infra 
red  Spectrograms. 
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the  case  in  the  present  work,  the  calculated  frequencies 
will  be  the  same  for  either  the  cis-or  trans-form  of  the 
molecule.    The  only  difference  will  be  in  the  selection 
rules  applying.    Thus  the  C2h  point  group  was  assumed 
as  applying  to  the  molecule  since  the  selection  rules  with 
this  assumption  permitted  better  agreement  between  the  more 
intense  observed  lines  and  calculated  lines.     The  matching 
required  between  the  calculated  accidental  degeneracies  and 
observed  lines  indicates  that  the  resonance  has  caused  a 
rather  lar^e  shift  in  the  position  of  these  observed  lines. 


TABLE  III 

Comparison  of  Calculated  and  Observed  Frequencies,  in  Wave 

Numbers 


Observed 
I.R. 


520 

909  v.v.S. 

1014-  v.s. 

1282  m. 

1374  S. 
1391  S. 

1479  w. 
1497  w. 
1597  v.v.S. 
1608  v.v.S. 

1759  w. 

1812  m. 
1821  m. 
1969  w. 

2495  w. 


Ra. 

320  v.w.O) 
513    S.  (P) 

686  v.w.(?) 

910  m(D) 

1205  vS(P) 
1279  vS(P) 
1302  S.(^) 


1442  v.S.(P) 

1643  v.v.S. (P) 
1780  w(P) 

2092  v.w.(?) 
2842  v.w.(?) 


Calculated 
434  Ra 

894  I.R. 

1180  Ra 

1230  I.R. 
1320  Ra 

1400  Ra 
1640  Ra 


1640  Ra 
1730  I.R. 


2062  Ra 
2430  I.R. 
2846  Ra 


TABLE  III  (Continued) 

Observed  Calculated 
2879  vn.w(?) 

2981  w.  2960  I.R. 

3014  v.S.(P)  2990  Ra 

^015  v.S.  3020  I.R. 
^095  v.S. 

3101  m(o)  2990  Ra 

4938  w  4930  I.R. 

The  form  of  the  normal  vibrations  can  be  de- 
termined4^ by  substituting  the  corresponding  value  of  \  in 
the  secular  equation  and  solving  for  the  ratios  of  the  co- 
ordinates which  are  equal  to  the  ratios  of  the  minors  of 
any  one  row  of  the  determinant.    These  ratios  always  remain 
equal  during  the  oscillation  of  the  particles  for  any  one 
normal  vibration.    For  A.^,  in  species  B^,  the  ratios  of  the 
coordinates 

yl  :  y2  5  y<*  :  y3  :  *1  :  s2  :  zoC  :  z3 
were  calculated  to  be  in  the  ratios 

.171063:  -.020901:  -.155850:     .950302:  -.812493 
:  1.03c836:  -.031805:  -.108816 


43.    Herzberg,  0d.  Git.,  p.  69 
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The  form  of  this  normal  vibration  is  shown  in  Fig.  3* 


Figure  3« 

The  solution  of  these  determinants  was  accom- 
plished by  Crout's  method4"4-  for  solving  determinants.  This 
method  is  also  readily  adaptable  to  computing  machines  and 
has  the  advantage  as  well  of  having  a  continual  check  on 
the  numerical  work.    The  forms  of  the  remaining  normal  co- 
ordinates could  be  obtained,  if  desired,  by  determining  the 
ratios  of  the  minors  of  any  row  of  both  of  the  determinants 
for  species  A-.  and  B-.  for  all  of  the  remaining  values  of  X. 


44.     Crout,  P.  D. .  Transactions  Am.   Inst.  Elec.  Eng.  IX  (1941) 


r 
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5.    An  Approximate  Calculation 
In  order  to  avoid  the  expansion,  of  the  eighth 
and  ninth  order    secular  equations  prior  to  the  adaption  of 
Bryan's  method,  an  approximation  scheme  was  used,  the  re- 
sults of  which  should  be  of  interest  in  a  comparison  with 
the  results  obtained  by  an  exact  determination.     The  ap- 

45 

proximation  scheme  that  was  used  was  suggested  by  Herzberg  ' , 
Consider  first  the  force  constants  associated 
with  the  bond  stretching  and  bending  of  the  hydrogen  atoms 
infinite  and  calculate  the  frequencies  for  the  carbon-hydro- 
gen groups  vibrating  as  one.     Then  consider  the  masses  of 
the  carbon  atoms  infinite  and  calculate  the  vibrations  of 
the  hydrogen  atoms  about  the  carbon  atoms. 

When  the  hydrogen  force  constants  are  consider- 
ed infinite,  the  only  symmetry  coordinates  that  remain  in  the 
potential  energy  expression  are  Ar^,  ^r23*  ^^34,  2HdA02« 
The  only  cartesian  coordinates  that  need  be  considered  then 
are  a^,  *2*  *3*  fr4*  ^2'  y3  and  y4*    For  sPecies  Ai 

these  are  related  by  the  symmetry  treatment  as  follows: 

M+2m  x 
*4      *!'  *3  =  *2  =  "  *1«  y4  =  "*H  y3=  "y2-  (1) 

This  then  requires  the  solution  of  only  a  third  order  deter- 
minant • 

Under  these  conditions  the  terms  in  the  poten- 
tial energy  expression  are 


45.    Herzberg,  Op.  Cit..  p.  195 


-  r      •  «c  r 


I^i2 12  =/|  (y2-y1-a)2+(-^vbJ2  -  V°i2\} 

|Ar23l2  =    [-2  y2+C  -|43  II  2 

1^12  -   i    f  p^-^x-2^ 


(2) 

(3) 

(4) 


G. 


2  _ 


I  (  UT(  y2-y i+a }  *+  ( *>iai+b }  *H ( - 2y2+ c  ^ 1 ' 

-cos©^   2  (5) 

i      ("  F  (-2y2+c)("yi+y2+a) 

sin^62o  (  J  [(-2y2+c)^]I  (-v^v^a^u^-b^] 


-cos© 


where 


=  1  - 


2  j 
M  -  2M 


(6) 


M  -  M 

The  expansion  of  the  secular  ecuation  gives 

a  third  order  equation  in  X.  which  is  then  solved  by  Horner's 

method  to  give  the  three  roots.     The  three  roots  for  this 

equation  give  three  normal  frequencies  of  vibration  which, 

expressed  in  wave  numbers,  are 
1 

H 

Now  to  obtain  the  frequencies  associated  with 
the  hydrogen  motion  the  carbon  masses  are  considered  infinite 
which  means  that  all  ccrrdinates  associated  with  motion  of 
the  carbon  atoms  disappear. 

The  only  remaining  coordinates  from  the  sym- 
metry conditions  are  then 


=  289cm'1  ,    ^  =  1140cm"1  ,  i    =  1450cm"1.  (7) 
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y5  =  -yr,  *6  -  *r  =  -  (**  +z3).  (8) 

This  leaves  five  coordinates  and  requires  the  solution  of  a 
fi^th  order  secular  equation.     The  five  roots  of  the  secular 
equation  again  give  five  normal  frequencies  of  vibration 
which  are 

J    =  882cm"1  ,  J    -  1350cm"1  ,  I    =  2840cm"1 
h  X2  X3 

=  2920cm'1  ,    I    =  3630cm"1.  (9) 

For  species        when  the  hydropen  bond  constants 
are  considered  infinite,  the  coordinates  remaining  from  the 
symmetry  treatment  are: 

y2  -  ~*Vi  and  z2  =  Ry"i  ~  Hzi 


where 


n  =     M+2m       R  -  X+  e  +»-n.       tt  =      f  - 1 

11       M+m    '  n  1+x      <  n  1+X 


and 


x  -  2(b+h)    p  -  2(b+e+h)     e  =  2a+2f+c     f_  2a-2d+c 


_  .     Mc  „      2M(b+h)          _  M(2a+c)  ,,nv 

^     m(2g+c;  *  Z1      m(2g+c)  -      "  m(2g+c)  *  llu; 

The  process  of  solving  for  the  roots  is  the 

same  as  that  outlined  for  species  A^.     In  this  case  the 

solution  of  a  second  order  secular  equation  is  required 

which  gives  the  frequencies 

£    =  1600cm"1      and    £    =  678cm"1.  (11) 

In  species  B-^  when  considering  the  frequencies 
due  to  the  motion  of  the  hydrogen  atoms  the  same  terms  dis- 
appear in  the  potential  function  and  the  kinetic  energy 
function  as  for       but  the  relations  between  the  remaining 


r  r 


cartesian  coordinates  are  as  follows: 

*is  =  $*^  -(fr'+9e)*"<t'*s**ir  -*#e+ey4- (12) 

This  then  requires  the  solution  of  a  fourth  order  secular 
equation  which  gives  the  frequencies 

I    =  1620cm"1,  J    ■  2050cm"1,  i    =  2770cm"1,  1    =  2830cm"1.  (13 

'    A-2  A,^  A.^ 

First  of  all,  it  is  noticed  in  Table  IV  that  in 
the  process  of  making  the  approximation  three  roots  dis- 
appear.   The  remaining  roots,  while  having  a  similar  overall 
pattern  to  those  determined  by  the  exact  solution,  differ 
from  them  by  about  200  wave  numbers  in  the  case  of  the 
carbon  oscillations  and  by  a  greater  amount  in  the  case  of 
the  hydrogen  oscillations.    The  wave  numbers  in  the  approxi- 
mation method  are  smaller  in  general  than  in  the  exact  de- 
termination as  would  be  expected  from  the  general  theory  of 
oscillations. 

TABLE  IV* 

Comparison  of  Frequencies,  in  Wave  Numbers 
"From  Approximate  and  Exact  Calculations 


Alc 

Alh 

Blc 

Bih 

434 

289 

... 

894 

678 

•  •  • 

1180 

... 

882 

1230 

•  •  • 

... 

1320 

1140 

... 

1400 

.  •  • 

1640 

1450 

... 

17^0 

1600 

. .  • 

1640 

1350 

2430 

... 

1620 

2060 

2960 

... 

2050 

2860 

2840 

3020 

... 

2770 

2990 

2920 

4930 

... 

2830 

2990 

^630 

(*  indicates  exact  calculations,  the  subscript  c  indicates 
carbon  vibrations,  and  h  indicates  hydrogen  vibrations) 


6 .    Calculation  of  Heat  Capacity 
From  a  knowledge  of  the  fundamental  frequencies 
of  a  molecule  it  is  possible  to  calculate  its  thermodynamic 
quantities.     It  is  generally  best  to  calculate  the  partition 
function  first  since  all  of  the  other  thermodynamic  quanti- 
ties can  be  determined  from  it.    The  partition  function  is 

,  46 
given  by 

(i) 

where  €n  is  the  total  energy  of  the  atoms  or  molecules  in  a 
given  equilibrium  distribution,  k  is  Boltzmann's  constant, 
T  is  the  absolute  temperature  and  gn  is  the  statistical 
weight  or  degeneracy.     This  can  be  separated  into  a  trans- 
lational  part  and  an  internal  part.     The  internal  partition 
function,  which  is  the  part  of  interest  here,  is  given  by 

QvrJL  '  JET  fur*,   e      ^  (2) 

This  part,  again,  is  composed  of  con- 
tributions from  electronic,  vibrational  and  rotational  energy, 
but  in  most  cases  the  electronic  contribution  is  negligible. 
Thus,  there  are  only  two  parts  remaining,  the  vibrational 
and  rotational  energy  contributions.    The  partition  function 
can  then  be  written  * 

where 

46.    Herzberg,  Op.  Cit..  p.  501  ff. 
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and  F  (J)  is  the  rotational  term  value,  J  is  the  rotational 
quantum  number,  GQ  is  the  vibrational  term  value,  and  v^  are 
the  vibrational  quantum  numbers.     If  the  interaction  between 
vibration  and  rotation  is  neglected,  this  can  be  written  as 
the  product  of  the  vibrational  and  rotational  partition 
functions,  and  each  Dart  can  be  separately  evaluated.  Then 

*Tow,  in  the  usual  harmonic  oscillator  aDproxi- 
mation  we  have 

Gr0  C^i  ,  ^  •  •  • )  =  21  oji  nri  (  ^ 

where  (J^  are  the  frequencies  in  cm""*"  for  infinitesimal 
amplitudes.     Then  h  ,  * 

I  at 

which  can  be  written  .       ,  , 

by  making  use  of  the-  relation 

If  degeneracy  exists  and  d^,  ^3***-  represent  the  degrees 

of  degeneracy  of  6J-.    <*-U   f»j    ...  the  partition  function  can 

be  written  _^  Jt^.-^.  ^       _      Xe-  N-^? 

^.fi-e  ^  ...  (8 


In  actual  calculations  the  wave  numbers        are  used  instead 
of  the  cu>^. 

For  the  rotational  partition  function  the  expres- 
sion varies  depending  upon  whether  the  molecule  is  linear, 
symmetric  top  or  asymmetric  top.    For  a  linear  molecule4"^ 
47.     Herzberg,  Op.  Cit..  p  505- 


I or  a  rigic  symmetrxc  luu  muiecuie 
and  for  asymmetric  top  molecules 


9) 

for  a  rigid  symmetric  top  molecule 

and  for  asymmetric  top  molecules 

where  A,  B  and  C  are  the  rotational  constants. 

Other  thermodynamic  quantities  can  then  be 
calculated  in  terms  of  the  partition  function.    For  example, 
the  specific  heat  at  constant  pressure  is  given  by 

Cp>uvck  =  njfLT  ■ — j^p    J  .  (12) 

Again,  if  the  interaction  of  rotation  and  vibration  is  neg- 
lected, then 

and 

<^r>nr     v-r        sit  (14) 

For  almost  all  gases  the  rotational  contribution  is  practi- 
cally equivalent  to  the  classical  values  for  ordinary  temper- 
atures.   By  substituting  the  value  for        In  the  above  equa- 
tion the  value  for  the  specific  heat  is  . 


(15) 

The  calculations  for  the  value  of  the  vibrational 
contributions  to  the  heat  capacity  are  simplified  by  the 


use  of  a  table  published  by  Wilson*  where  the  vibrational 
contributions  to  the  heat  capacity  are  given  as  a  function 
of  •  The  calculated  value  for  Butadiene  determined  for 
a  temperature  of  25°C  with  the  aid  of  this  table  is 

7.26  cal/mol  deg. 
By  adding  this  to  the  translational  and  rotational  contri- 
butions as  given  by  Wilson  the  heat  capacity  turns  out  to 

be  _q1 

1^.21  cai/mol  deg. 

The  best  experimental  value  available  obtained  by  comparison 

of  the  results  published  by  Astor,  Moessan,  Gustave,  Hardy, 

4.0  crn 
and  Szasz      and  Scott  and  Mellors  is 

19.16  cal/niol  deg. 
The  difference  between  the  observed  value  and  the  calculated 
value  is  due  to  internal  rotation  and  could  aid  in  estab- 
lishing the  height  of  the  potential  barrier  in  the  molecule. 
7.    Evaluation  of  Results 

The  results  obtained  in  this  calculation  indi- 
cate again  the  fact  that  it  is  oossible  to  transfer  force 
constants  from  parts  of  one  molecule  to  another  where  similar 
structure  exists.     The  fairly  reasonable  check  of  the  pat- 
tern of  the  observed  spectra  with  the  calculated  lines  de- 
monstrates that  the  force  constants  as  chosen  in  the  poten- 
tial function  probably  agree  favorably  with  the  correct 

48.  Wilson,  E.  B.  Jr.,  Chem.  Rev.  27,  17  (1940) 

49.  Astor,  John  G.  Moessan,  Gustave  W. ,  Hardy,  Howard  C, 
and  Szasz,  George  J.     J.  Chem.  Phys.  12,  458  (1944) 

50.  Scott,  Russell  B.  and  Mellors,  Jane  W. ,  J.  Res.  Nat. 
Bu.  Stand.  34,  243  (1945) 
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values.    The  characteristic  frequencies  of  various  groups 
within  molecules  and  of  the  motion  within  these  groups  as 
given  by  Herzberg^    check  fairly  well  with  calculated  fre- 
quencies which  serves  as  an  aid  in  identifying  certain  mo- 
tions in  the  molecule.    Exact  agreement  between  observed 
and  calculated  lines  is,  of  course,  not  possible  because 
the  valence  force  model  is  an  approximation,  and  there  may 
also  be  complications  due  to  impurities  in  the  experimental 
samples  and  because  of  the  prescence  of  combination  and  over 
tone  bands. 

Probably  of  greater  significance  in  this  work  is 
the  method  used  for  expanding  the  secular  equation.  This 
part  of  the  problem  has  always  been  the  most  tedious  and 
laborious  of  all  the  work  involved.     By  the  use  of  this 
method,  not  only  is  it  uneccessary  to  make  use  of  approxi- 
mations, but  the  exact  expansion  can  in  general  be  accomp- 
lished in  less  time  and  with  less  effort  than  previous  ap- 
proximation schemes  required.    Furthermore,  a  continuous 
check  is  provided  on  the  numerical  work  thereby  reducing 
the  chance  of  error.    Since  the  method  is  also  readily  adap- 
table to  I.B.M.  calculation,  an  additional  saving  in  time 
and  labor  can  be  accomplished  without  the  use  of  iteration 
or  approximation  schemes  with  their  attendent  loss  or  shift 
of  frequencies. 

51.      Herzberg,  0p_.  Cit . ,  p.  195 
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ABSTRACT 

The  study  of  infra-red  and  Raman  spectra  and  their 
relationship  to  the  vibrations  and  rotations  of  atoms  in 
molecules  has  become  more  and  more  a  matter  of  interest  to 
the  physicist  and  chemist  within  the  past  several  decades. 
Thp  study  of  the  vibrations  of  the  atoms  which  give  rise  to 
the  broad  band  spectra  in  the  near  infra-red  has  led  tc  a 
more  thorough  understanding  of  the  forces  in  the  molecules 
including  a  theoretical  explanation  of  valence.    The  study 
of  the  rotation  of  groups  of  atoms  within  a  molecule  as  well 
as  the  rotation  of  molecules  as  a  whole  which  give  rise  to 
the  band  spectra  in  the  far  infra-red  and  to  fine  structure 
in  the  near  infra-red  has  led  to  a  rather  accurate  means  for 
determining  the  structure  of  molecules. 

This  paper  deals  with  the  theoretical  determination  of 
the  vibrational  frequencies  of  1,  3-  butadiene.    As  is  well 
known,  the  frequencies  obtained  from  a  quantum  mechanical 
treatment  of  the  problem  are  the  same  as  those  resulting  from 
a  classical  determination;  consequently,  though  the  quantum 
mechanical  methods  are  outlined,  the  actual  calculations  are 
carried  out  by  the  classical  method  of  treating  small  oscil- 
lations.    The  normal  frequencies  of  vibration  determined  for 
an  oscillating  system  having  particles  of  masses  equal  to 
those  of  the  masses  of  the  atoms  with  force  constants  of  mag- 
nitude equal  to  the  forces  actually  existing  in  a  molecule 
will  be  the  same  as  the  infra-red  and  Raman  frequencies 
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emitted    or  absorbed  by  the  molecule. 

In  general  the  force  constants  are  unknown  so  that  the 
problem  essentially  has  to  be  worked  backward  by  carrying  it 
out  algebraically  and  then,  by  making  use  of  obsprved  spectra, 
determining  the  force  constants.    For  a  large  molecule,  such 
as  butadiene,  this  is  not  easy  nor  even  feasible.  Consequently, 
force  constants  determined  by  others  for  smaller  molecules 
where  the  structure  of  part  of  the  molecule  is  similar  to  a 
part  of  butadiene  are  used  for  a  determination  of  the  normal 
frequencies.     The  force  constants  in  this  work  were  transfer- 
red from  the  work  on  propylene  by  Wilson  and  Wells''"  and  from 

a  paper  on  force  constants  in  organic  molecules  by  Crawford 
p 

and  Brinkley  . 

The  cis-form  was  selected  as  the  normal  configuration 
for  butadiene.    The  force  field  chosen  was  the  valence  type 
field.    This  means  that  the  changes  in  bond  length  and  changes 
in  the  bond  angles  from  the  equilibrium  configuration  were 
chosen  as  the  coordinates  for  the  potential  function.  The 
application  of  group  theory  to  the  symmetry  in  the  molecule 
factored  the  secular  equation  obtained  ^rom  the  coefficients 
of  the  variables  in  the  kinetic  and  potential  energy  expres- 
sions from  a  determinant  of  order  twenty  four  to  four  deter- 
minants of  order  nine,  eight,  four  and  three.     Since  the 

1.  Wilson,  E.B.  Jr.,  and  Wells,  A.J.,  J.  Chem.  Phy.  9,  314 
(1941) 

2.  Crawford,  B.L.Jr.,  and  Brinkley,  S.R.  Jr.,  J.  Chem.  Phy. 
9,  69  (1941 
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latter  two  contained  only  out-of-plane  vibrations  and  would 
thus  give  zero  frequencies  with  the  type  of  **orce  field 
chosen,  only  the  two  equations  of  order  nine  and  eight  had 
to  be  solved.     These  are,  however,  in  determinantal  form, 
and  the  process  of  expanding  secular  equations  of  this  size 
is,  in  general,  very  laborious  and  tedious. 

In  order  to  reduce  the  amount  of  labor  involved,  an 
approximation  scheme  was  attempted.    As  a  first  approxi- 
mation, the  carbon-hydrogen  bond  force  constants,  both 
stretching  and  bending,  were  considered  infinite.  This 
served  then  to  increase  the  mass  of  thp  carbon  atoms 
slightly  while  eliminating  the  hydrogen  frequencies.  Two 
secular  equations  of  order  three  and  two  had  to  be  expanded, 
which  yielded  an  approximate  value  for  the  carbon  skeletal 
vibrational  frequencies.     Secondly,  the  carbon  masses  were 
considered  infinite  and  then  the  motions  of  the  hydrogen 
atoms  with  respect  to  the  carbon  atoms  were  determined.  This 
required  the  expansion  of  two  secular  equations  of  order 
five  and  four. 

In  the  process  of  making  these  approximations,  three 
frequencies  disappeared,  which  should  have  or  would  have 
bpen  determined  by  an  exac+  expansion  of  the  original  secu- 
lar equations.     A  method  was  then  used  for  an  exact  expan- 
sion  of  the  secular  equations  which  was  due  to  J.  G.  Bryan^ 
and  was  developed  for  expanding  matrix  characteristic 

3.    Bryan,  J.G.,  A  Method  for  the  Exact  Determination  o^  the 
Characteristic  Equation  and  Latent  Vectors  of  a  Matrix 
with  Applications  to  the  Discriminant  ^unction  for  More 
than  Tv^o  Groups.  Harvard  Graduate  School  of  Education 
Dissertation  (1950) 
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equations  and  determining  the  latent  vectors  or  character- 
istic vectors  of  a  characteristic  equation.     The  character- 
istic equation  of  a  matrix  is  defined  by 

|A  *  \I|  -  0 

where  A  is  the  original  matrix  and  I  is  the  unit  matrix. 
While  the  secular  equation  is  not,  in  general,  in  the  form 
defined  above,  it  can  be  put  into  this  ^orm  by  adding  and 
subtracting  rows  and  columns  or  multiples  of  rows  and 
columns  until  the  variable  terms  appear  only  on  fhe  diagonal 
and  then  dividing  each  column  by  the  coefficient  of  the  vari- 
able in  that  column.     Since  the  secular  equation  will  then 
be  in  the  form  as  defined  above,  it  is  possible  to  make  use 
of  Bryan's  method  ^or  an  exact  expansion  of  the  characteristic 
equation.     The  method  is  described  completely  and  the  expan- 
sion of  the  eighth  order  secular  equation  is  carried  out  in 
detail  as  an  illustration  of  the  method.    The  actual  expansion 
of  the  secular  equations  by  this  process  requires  less  time 
and  effort,  in  general,  than  the  approximation  scheme.  This 
method  has  an  additional  advantage  in  that  a  continual  check 
is  provided  on  the  numerical  work  as  the  expansion  proceeds. 
Furthermore,  the  method  is  readily  adaptable  to  I.B.M.  auto- 
matic punch  card  computors  which  can  result  in  an  additional 
saving  in  time  and  effort  while  at  the  same  time  carrying  out 
an  exact  expansion  rather  than  an  approximation  with  its  at- 
tendant loss  or  shift  of  frequencies. 

After  the  equations  were  expanded,  they  were  solved  by 
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Horner's  method  yielding  the  seventeen  required  roots. 
These  roots  were  proportional  to  the  normal  frequencies  of 
vibration  of  the  molecule  and  after  solving  for  the  fre- 
quencies in  terms  of  wave  numbers,  they  were  compared  with 
observed  frequencies.    The  comparison  showed  a  fairly  good 
agreement  between  the  observed  and  calculated  frequencies. 
A  comparison  was  also  made  between  the  frequencies  obtained 
by  the  exact  expansion  and  the  approximation.    This  showed 
that  the  carbon  vibrations  in  the  approximation  were  all 
fairly  uniformly  displaced  by  about  two  hundred  wave  numbers 
toward  lower  frequencies  as  is  to  be  expected  from  the  gener- 
al theory  of  small  oscillations.     The  hydrogen  frequencies  were 
displaced  by  an  even  larger  amount  as  also  was  to  be  expected 
but  not  in  as  regular  a  manner  as  the  carbon  oscillations. 

The  results  of  this  paper  again  indicate  that  it  is 
possible  to  transfer  force  constants  from  parts  of  one  mole- 
cule to  another  where  they  exist  in  similar  surroundings. 
The  fairly  good  correlation  between  the  observed  and  calcu- 
lated frequencies  indicates  that  the  force  constants  as  chosen 
are  very  nearly  the  correct  values. 

Probably  of  more  significance  in  this  paper  is  the  appli- 
cation of  Bryan's  method  to  the  exact  expansion  of  the  sec- 
ular equation.    This  eliminates  the  necessity  of  introducing 
approximations,  and  the  expansion  can  be  carried  out  with  less 
work  and  more  of  a  check  on  the  work  than  previously.  This 
expansion  as  well  as  the  solution  of  the  expanded  equation  can 
be  carried  out  on  I.B.M.  computors.    Furthermore,  the  deter- 


mination  of  the  normal  modes,  which  requires  the  solution  of 

the  ratios  of  the  minors  of  any  row  of  the  secular  equation 

with  a  given  value  for  the  frequency  substituted,  was  per- 

k 

formed  by  Crout's    method  for  the  solution  of  determinants 
which  is  also  readily  adaptable  to  I.B.M.  computors. 


l+o      Crout,  P.D.,  Trans.  Am.  Inst.  Elec.  Eng. ,  LX  (19kl) 
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